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I.     INTRODUCTION 

The  Naval  Undersea  Warfare  Engineering  Station  (NUWES)  conducts  testing  of 
torpedoes  on  various  underwater  ranges  in  the  Pacific  Northwest  (Washington  State  and 
Canada).  On  occasion  these  torpedoes  malfunction  and  are  lost  on  or  under  the  ocean 
floor  and  must  be  recovered.  The  recovery  operation  involves  locating  the  torpedo,  if 
neccessary  digging  it  out  from  under  the  ocean  bottom,  and  attaching  a  cable  to  it  in 
order  to  pull  it  out.  This  is  all  done  by  remote  control  by  an  operator  on  a  surface  ship 
using  video  and  sonar  images  to  guide  him.  The  sonar  signal  is  used  for  general  position 
location  and  then  video  images  are  used  to  conduct  the  remainder  of  the  operation. 

There  are  several  problems  which  occur  during  the  acquisition  of  these  images. 
Probably  the  most  critical  of  these  are  caused  by  the  artificial  lighting  required  at  depths 
of  up  to  1500  feet.  Uneven  lighting  can  produce  severe  shadows  and  bright  blooming 
on  the  same  image  frame.  Bottom  texture  (important  for  finding  the  entry  point  when 
the  torpedo  has  buried  itself)  often  cannot  be  distinguished.  Silt  and  debris  thrown  up 
by  the  digging  device  can  totally  obscure  vision  for  up  to  several  hours  at  a  time. 

The  tools  chosen  to  help  solve  these  problems  were  an  IBM  AT  Personal 
Computed  with  a  PCVlSIONplus  FRAMEGRABBER2  board  installed  and  the  ITEX 
PCplus^  library  of  image  processing  subroutines.  This  system  allows  one  to  build  and 
run  software  tailored  to  his  her  particular  image  processing  needs.  The  system  is  rela- 
tively inexpensive  (when  compared  to  more  sophisticated  image  processing  work- 
stations), portable,  and   easily  installed  in  an  office  or  on  board  a  recovery  vessel. 

The  purpose  of  this  thesis  was  to  develop  an  image  processing  package,  using  the 
tools  described  above,  to  provide  contrast  and  other  enhancement  of  video  images.  Be- 
cause the  images  to  be  processed  are  used  as  feedback  in  operating  a  remote  control 
device,  a  real-time  or  near  real-time  approach  to  the  problem  was  neccesitated.  In  a 
previous  thesis,  [Ref.  1]  Roberto  Ventura  suggested  the  use  of  lookup  table  modification 
as  a  primary  solution.  There  were  also  several  other  non-lookup  table  algorithms  which 
appeared  to  have  potential  usefulness.    Once  developed,  the  package  was  to  be  tested 


1  Trademark  of  IBM 

2  Trademark  of  Imaging  Technology  Inc. 

3  Trademark  of  Immune  Technolosv  Inc. 


at  NUWES  and  delivered  to  them  for  their  use.  Inherent  in  this  requirement  was  a 
"user-friendly"  presentation  at  the  computer  monitor  and  appropriate  documentation  for 
its  operation.  The  text  of  this  thesis  shows  an  analysis  of  the  problem,  describes,  in 
general,  the  tools  used,  and  outlines  the  development  of  the  final  package.  Further,  it 
describes  the  results  of  the  testing  conducted  on  a  recovery  vessel  at  NUWES.  The  ap- 
pendices include  instructions  for  the  use  of  the  package  and  the  C  language  code  for  the 
program. 


II.     ANALYSIS  OF  THE  PROBLEM 

After  observing  a  recovery  operation,  receiving  comments  from  operators,  and  ex- 
amining video  tape  of  typical  underwater  conditions,  three  viewing  regimes  were  identi- 
fied. These  were  viewing  through  turbid  water  (due  to  disturbance  of  bottom  silt),  clear 
water  viewing  of  the  bottom  while  searching,  and  clear  water  viewing  of  shiny  objects 
which  may  cause  blooming.  Each  regime  was  analyzed  to  determine  what  problems,  if 
any,  could  potentially  be  solved.  The  calculation  and  display  of  histograms,  a  primary- 
tool  in  the  analysis,  will  be  discussed  in  a  later  chapter. 

The  most  challenging  regime  is  that  of  viewing  objects  through  turbid  water.  Often 
the  density  of  particles  in  the  water  is  so  great  that  all  incident  light  is  scattered  or  re- 
flected before  it  can  reach  the  area  or  object  to  be  observed.  In  this  case  the  video  sys- 
tem is  of  no  use  and  some  other  portion  of  the  electromagnetic  spectrum  should  be 
considered  for  image  acquisition.  However,  if  even  a  little  light  is  reflected  from  the 
object  of  interest  back  to  the  camera  there  is  potential  for  enhancement  of  the  resulting 
image.  The  histogram  of  a  typical  turbid  water  image,  as  shown  in  Figure  1,  provides 
information  which  could  lead  to  the  method  used  to  acheive  this  improvement.  The 
tight  distribution  of  pixel  values  around  the  center  range  of  intensity  levels  readily  lends 
itself  to  contrast  enhancement  using  look-up  table  modification  techniques. 

When  searching  for  a  torpedo  that  has  buried  itself  under  the  bottom,  the  general 
location  is  first  determined  using  sonar.  Once  in  the  vicinity,  however,  it  is  very  useful 
if  the  entry  point  can  be  found.  This  is  often  difficult  since  the  mud  closes  in  over  the 
top  of  the  torpedo.  It  would  be  helpful  if  some  method  could  be  found  to  better  identify 
marks  on  the  bottom  which  could  be  entry  points.  Figure  2  shows  the  histogram  of  a 
typical  image  of  the  bottom.  Again,  most  pixel  values  are  concentrated  around  the 
center  range  of  intensity  levels  and  this  provides  excellent  opportunity  for  enhancement 
by  look-up  table  modification.  Also,  an  enlarged  view  of  a  portion  of  the  image  could 
be  potentially  helpful  in  object  identification.  This  could  aid  in  identifying  a  partially 
uncovered  object  during  the  digging  operation. 

Blooming  is  a  phenomona  caused  by  the  specular  reflection  of  bright  lights  from 
shiny  objects.  It  often  occurs  after  the  torpedo  is  found  while  attempting  to  attach  a 
cable  to  it.  The  result  is  a  "white  out"  over  a  significant  portion  of  the  image.  The 
current  operational  solution  to  this  problem  is  to  dim  the  entire  image  on  the  monitor 


Figure  1.      Histogram  of  a  Turbid  Water  Image 


Figure  2.      Histogram  of  a  Bottom  Image 


using  the  brightness  control.  This  reduces  the  detail  available  from  the  darker  portions 
of  the  image.  The  histogram  of  an  image  with  blooming,  shown  in  Figure  3,  also  sug- 
gests that  a  look-up  table  modification  could  be  performed  to  reduce  the  brightness  of 
higher  intensity  pixels  without  dimming  the  lower  values. 


Figure  3.      Histogram  of  Image  With  Blooming. 


III.     DESCRIPTION  OF  THE  IMAGE  PROCESSING  SYSTEM 

A.  GENERAL 

This  chapter  describes  the  tools  to  be  used  in  the  attempt  to  improve  the  quality  of 
the  underwater  images.  The  image  processing  system,  shown  in  Figure  4,  includes 
hardware  in  the  form  of  a  personal  computer  with  the  I  TEX  PCplus 
FRAMEGRABBER  installed,  and  a  library  of  image  processing  software  called  ITEX 
PCplus. 

B.  HARDWARE 

The  framegrabber  (enclosed  within  the  dashed  line  in  Figure  4)  is  the  heart  of  the 
image  processing  system.  It  is  installed  in  a  host  computer  for  control  of  its  components 
and  for  disk  access.  Compatibility  with  the  host  is  ensured  by  placement  of  jumpers  in 
the  proper  locations  on  the  board.  This  also  sets  the  memory  base  and  register  base 
addresses.  The  framegrabber  receives  video  input  from  a  source  (camera  or  video  tape 
player)  and  sends  video  output  to  a  monitor.  At  the  input  the  image  is  digitized  by  an 
analog-to-digital  converter.  It  is  then  sent  through  the  input  lookup  table  (LUT).  There 
are  eight  hardware  lookup  tables  available  on  the  input  side.  One  is  selected  for  use  at 
any  given  time  and  can  be  loaded  with  the  desired  values  from  the  computer.  Values 
may  also  be  read  from  the  selected  lookup  table  and  stored  to  disk.  After  passing 
through  the  input  lookup  table  the  image  is  placed  into  one  of  two  memories,  A  or  B, 
whichever  is  selected  from  the  computer.  The  contents  of  these  memories  are  accessible 
from  the  computer  for  analysis,  processing,  or  storage  to  disk.  Once  an  image  is  in  a 
selected  memory,  it  may  be  sent  for  display.  Enroute  it  passes  through  one  of  eight 
output  lookup  tables  which  operate  in  the  same  manner  as  the  input  lookup  tables. 
Finally  it  is  converted  to  an  analog  image  to  be  displayed  on  the  video  monitor.  In  the 
case  of  a  continuous  video  image  from  the  source  this  entire  process  may  occur  even' 
one-thirtieth  of  a  second,  presenting  a  continuous  image  to  be  displayed  at  the  output. 
Single  frames  from  a  source  or  from  the  computer  disk  may  also  be  placed  into  memory 
and  presented  at  the  output  for  display.   [Ref.  2] 

C.  SOFTWARE 

The  framegrabber  hardware  is  controlled  using  a  library  of  functions  called  ITEX 
PCplus.    These  functions  may  be  incorporated  in  a  program  using  the  "C"  language. 
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Figure  4.      The  Image  Processing  System 


They  perform  various  tasks  including  computer  setup,  image  acquisition,  image  storage 
to  disk,  lookup  table  operations,  and  single  frame  image  processing.  Along  with  basic 
functions  to  place  the  framegrabber  into  operation,  several  others  were  identified  as 
having  potential  for  providing  real  time  or  near  real  time  contrast  enhancement.  These 
building  blocks  are  described  here. 

The  setup  functions  configure  the  computer  to  be  compatible  with  the  framegrabber 
configuration,  select  memory  size  options,  initialize  the  framegrabber,  and  place  it  into 
operation.  They  also  allow  selection  of  memory  for  acquisition,  display,  storage,  and 
processing.   A  brief  description  of  these  functions  follows: 

sethdw  Defines  the  hardware  settings  including  register  base  address,  memory  base 
address,  and  memory  option.  Register  base  and  memory  base  address  ar- 
guments tell  the  computer  how  the  framegrabber  is  configured.  The  mem- 
ory option  defines  either  dual  512x512  frame  memories  or  a  single  1024x512 
memory. 

setdim  Defines  the  memory  size.  Sets  memory's  horizontal  and  vertical  size  and 
number  of  bits  used  to  store  each  pixel. 


initialize  Initializes  all  hardware.  Sets  all  hardware  registers  and  lookup  tables  to 
standard  values. 

fgon  Turns  on  the  framegrabber  memory. 

select_mem  Selects  frame  memory.  Allows  selection  of  memory  for  loading  from  disk, 
storage  to  disk,  or  processing.  May  be  used  to  access  one  memory  while 
another  memory  is  being  displayed. 

display jmem  Selects  frame  memory  display.  Allows  selection  of  memory  for  acquisi- 
tion from  source  and  or  display  to  screen. 

fgoff  Turns  off  the  framegrabber  memory.  [Ref.  3] 

The  image  acquisition  functions  allow  input  of  image  data  from  the  source.  They 
allow  continuous  passage  of  image  frames  from  source  to  video  monitor,  single  frame 
storage  and  display,  and  memory  clear  operations.   They  are  described  as  follows: 

snap  Acquires  a  single  image  frame  from  the  source,  stores  it  in  frame  memory, 

and  displays  it  on  the  monitor.  Allows  single  frame  processing  or  storage 
to  disk  of  an  image  from  the  source. 

grab  Acquires  and  displays  image  frames  continuously  up  to  a  specified  number 

of  frames  or  indefinitely  until  new  a  command  is  given.  Allows  real  time 
imaee  processing  since  image  frames  pass  continuouslv  through  lookup  ta- 
blesT 

sclear  Clears  entire  screen  to  a  value.  Erases  frame  memory  and  screen  by  placing 
a  specified  intensity  level  into  each  memory  location. 

saveim  Saves  a  specified  portion  of  an  image  to  a  file  on  disk.  The  file  name  is 
passed  as  an  argument  to  the  function. 

readim  Reads  an  image  file  from  disk.  The  file  name  is  passed  as  an  argument  to 
the  function.   [Ref.  3] 

The  lookup  table  functions  provide  the  capability  to  select  a  specified  lookup  table 
for  use,  read  lookup  table  values  from  disk,  write  lookup  table  values  to  disk,  and  modify 
lookup  tables  in  various  ways.   These  functions  are  described  below: 

setlut  Selects  which  of  eight  lookup  tables  on  the  input  or  output  bank  will  be 
active  and  available  for  storage  of  values. 

walut         Writes  an  array  of  values  to  a  selected  lookup  table. 

ralut  Reads  selected  lookup  table  values  into  an  array. 

savejut     Saves  the  values  in  a  selected  lookup  table  to  disk. 

read_lut  Reads  lookup  table  values  from  disk  and  places  them  into  the  lookup  table 
from  which  they  were  stored. 

linlut         Places  linear  lookup  table  values  into  selected  lookup  table. 

imlut         Places  inverse  linear  lookup  table  values  into  selected  lookup  table. 


loglut         Places  logarithmic  lookup  table  values  into  selected  lookup  table. 

abslut  Places  values  into  selected  lookup  table  which  invert  all  inputs  greater  than 
127. 

histeq  Performs  histogram  equalization  on  a  specified  portion  of  an  image  and 
places  computed  values  into  selected  lookup  table.  [Ref.  3] 

Several  non-real  time  image  processing  functions  were  also  considered.     These 

functions  work  on  images  fixed  in  frame  memory  and  perform  filtering,  pixel  replication, 

averaging,  and  histogram  operations.   They  are  described  here: 

sobel  Performs  a  Sobel  edge  detector  operation  on  a  specified  portion  of  an  im- 

age.   Displays  prominent  edges  of  the  image. 

sharpen     Sharpens  a  specified  portion  of  an  image  by  selectively  amplifying  edges. 

lopass  Performs  a  low-pass  filtering  operation  on  a  specified  portion  of  an  image 
and  displays  result. 

zoom  Enlarges  by  a  factor  of  two  the  portion  of  the  image  in  the  upper  left 

quarter  of  the  display. 

average  Acquires  the  selected  portion  of  a  specified  number  of  live  image  frames, 
stores  them  on  disk,  computes  the  average,  and  displays  results. 

histogram  Computes  a  histogram  of  a  specified  portion  of  an  image  and  places  values 
into  an  array.  Allows  selection  of  spatial  sampling  rate  in  vertical  and 
horizontal  directions  and  size  of  storage  bins  (array). 

pan  Shifts  the  entire  image  horizontally  by  a  specified  amount. 

scroll         Shifts  the  entire  image  vertically  by  a  specified  amount. 

rectangle  Draws  a  rectangle  with  specified  size  and  location  on  the  screen  and  in 
frame  memory.   May  be  superimposed  on  a  fixed  image.   [Ref.  3] 


IV.     DEVELOPMENT  OF  THE  MENU  PROGRAM 

A.  GENERAL 

The  concept  of  the  "Menu"  program  was  to  place  the  various  image  processing 
options  and  the  analytical  tools  into  a  format  which  allows  option  and  parameter 
changes  to  be  easily  made  and  evaluated.  Further,  the  menu  format  is  used  in  the  final 
image  processing  package  to  be  delivered  to  NUWES. 

This  chapter  describes  how  the  various  ITEX  PCplus  functions  were  incorporated 
using  Microsoft*  C  5.0  [Ref.  4]  to  produce  the  "Menu"  program.  It  outlines  how  the 
setup  functions  were  used  to  bring  the  board  on  line  and  perform  some  basic  operations. 
Then  it  describes  how  the  histogram  and  lookup  table  read  write  functions  were  used  to 
build  the  analytical  tools.  Finally,  it  shows  how  various  lookup  table  modification  and 
other  image  processing  functions  were  written  into  the  program  and  tested  for  applica- 
bility to  the  stated  problems. 

B.  BASIC  OPERATIONS 

1.  Setup 

Four  functions  are  required  to  define  and  initialize  the  hardware  and  place  the 
board  into  operation.  The  sethdw  function  is  used  with  arguments  which  define  the 
memory  base  address  as  D0000.  the  register  base  address  as  100,  and  select  the  dual 
memory  model.  It  is  also  used  in  a  menu  option  which  allows  the  operator  to  interac- 
tively enter  different  addresses  if  the  jumpers  on  the  framegrabber  board  are  configured 
differently.  The  setdim  function  is  used  to  set  the  frame  memory  size  to  512  x  512  pixel 
locations  with  8  bits  for  each  pixel  value.  The  initialize  function  is  included  in  the  pro- 
gram to  set  all  registers  and  lookup  tables  to  standard  values  and  fgon  is  included  to  turn 
the  framegrabber  memory  on. 

2.  Image  Acquisition 

Several  menu  options  were  devoted  to  image  aquisition  functions.  Two  options 
use  the  display_mem  function  to  allow  the  operator  to  choose  which  frame  memory  to 
display.  Two  other  options  use  the  select_mem  function  to  allow  selection  of  a  frame 
memory  for  processing  or  storage  to  disk.  Combinations  of  these  four  choices  permit 
one  image  to  be  displayed  while  another  is  being  processed.   The  snap  function  was  used 


4  Trademark  of  Microsoft  Corporation 


10 


in  a  menu  option  to  allow  the  acquisition  of  a  single  image  frame  into  a  selected  mem- 
ory. Another  option  selects  continuous  acquisition  and  display  of  image  frames  using 
the  grab  function.  The  clear  screen  menu  options  use  the  sclear  function  to  erase  the 
screen  and  frame  memory  to  either  black  or  white.  Finally,  access  to  images  stored  on 
disk  is  accomplished  with  menu  options  using  the  saveim  and  readim  functions. 

C.     BUILDING  THE  ANALYSIS  TOOLS 

1.     Computation  and  Display  of  the  Histogram 

The  histogram  was  the  most  important  tool  used  in  analysing  the  characteristics 
of  an  image.  The  histogram  function  provides  the  computed  histogram  values  in  an  ar- 
ray, but  methods  needed  to  be  devised  to  specify  the  portion  of  the  image  to  include  in 
the  calculation,  to  display  the  histogram  graphically,  and  to  compare  images  and  their 
histograms  before  and  after  processing. 

To  specify  the  portion  of  the  image  to  be  used  in  the  histogram  calculation,  a 
rectangle  whose  size  and  position  could  be  changed  from  the  keyboard  was  superim- 
posed on  a  fixed  image  displayed  on  the  monitor  using  the  rectangle  function.  This  se- 
quence (the  snap  and  rectangle  functions)  was  placed  inside  a  loop  to  give  the 
appearance  of  continuous  acquisition  and  display.  The  size  and  location  parameters  of 
the  rectangle  were  then  used  as  the  area  arguments  for  the  histogram  function.  A  default 
option  to  use  the  entire  screen  was  also  programmed  so  that  the  rectangle  routine  could 
be  bypassed  if  desired. 

A  method  was  also  needed  to  display  the  histogram  values  graphically  on  the 
computer  monitor  at  the  same  time  the  image  was  being  displayed  on  the  video  monitor. 
This  was  accomplished  in  a  routine  which  used  Microsoft  C  5.0  linedrawing  and  other 
graphics  functions  to  draw  vertical  lines  from  a  horizontal  axis  up  to  scaled  histogram 
values.  The  result  was  a  256  line  wide  (one  line  for  each  intensity  level)  bar  graph 
showing  the  relative  density  of  each  intensity  level  of  the  image.  Examples  of  the 
histogram  display  are  shown  in  Chapter  2  (Figures  1,  2,  and  3). 

To  make  a  comparison  of  an  image  and  its  histogram  before  and  after  process- 
ing requires  the  use  of  both  frame  memories.  An  image  frame  is  acquired  into  one 
memory  using  the  snap  function  with  a  linear  input  lookup  table.  The  other  memory 
is  selected  (the  display_mem  function)  and  new  values  are  loaded  into  the  input  lookup 
table  using  the  walut  function.  Then  another  image  frame  is  acquired.  The  time  differ- 
ence between  the  two  frames  is  so  short  that  the  two  images  are  virtually  the  same.  The 
histogram  values  for  each  frame  are  then  computed  with  the  histogram  function  and 
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stored  in  two  arrays.  The  histogram  and  frame  memory  to  be  displayed  can  then  be 
selected  from  the  keyboard,  presenting  processed  and  unprocessed  versions  of  the  image 
and  histogram  to  be  examined. 

2.     Manual  Modification  and  Display  of  Lookup  Tables 

As  stated  earlier,  the  lookup  table  is  the  primary  means  by  which  the  contrast 
of  the  image  frames  can  be  enhanced  in  real  time.  In  order  to  predict  and  evaluate  the 
effect  of  a  particular  set  of  lookup  table  values  on  an  image  a  means  of  displaying  those 
values  graphically  was  needed.  Further,  it  was  desired  to  be  able  to  change  the  lookup 
table  graph  interactively  from  the  computer  keyboard  and  evaluate  the  effect  of  these 
changes  on  the  image  frames. 

To  plot  the  lookup  table  the  ralut  function  is  used  to  read  the  values  from  the 
selected  lookup  table  into  an  array.  Microsoft  C  5.0  linedrawing  functions  with  line 
segment  endpoints  defined  by  array  values  plot  a  continuous  line  which  represents  the 
lookup  table.    Figure  5  shows  a  typical  lookup  table  display. 

Two  methods  were  programmed  to  interactively  modify  lookup  tables.  The  first 
method  reverses  the  process  used  to  plot  the  lookup  table  described  above.  Line  seg- 
ment end  points  are  entered  from  the  computer  keyboard.  The  line  is  then  drawn  on  the 
computer  monitor.  The  equation  for  that  line  is  computed  and  values  are  derived  to  be 
placed  in  the  appropriate  positions  in  the  lookup  table  array.  The  walut  function  then 
writes  the  array  to  the  lookup  table.  The  effect  on  the  image  is  immediately  evident  and 
may  be  evaluated.  The  second  modification  routine  interactively  varies  the  slope  of  the 
lookup  table.  The  point  where  the  line  intercepts  the  top  of  the  graph  is  shifted  to  the 
left  or  right  by  pressing  the  arrow  keys  on  the  computer  keyboard.  The  remaining  points 
on  the  line  to  the  bottom  of  the  graph  are  then  shifted  in  the  same  direction,  with  the 
amount  of  shift  varying  linearly  to  zero  from  the  initial  shift.  The  shape  of  the  line  re- 
mains very  close  to  the  original,  but  the  overall  slope  is  changed,  as  shown  in  Figures  6 
and  7.  As  the  point  by  point  changes  are  made,  the  lookup  table  array  is  continually 
updated  and  the  values  written  to  the  selected  lookup  table.  The  effect  of  each  change 
on  the  incoming  image  frames  occurs  immediately.  The  slope  may  also  be  varied  in  a 
similar  manner  by  shifting  the  bottom  portion  of  the  lookup  table  line. 

At  this  point  the  savejut  and  read_lut  functions  were  included  as  menu  options 
to  allow  future  use  of  particularly  effective  sets  of  lookup  table  derived  during  exper- 
imentation. 
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Figure  5.      Typical  Lookup  Table  Display 
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Figure  6.      Lookup  Table  Before  Slope  Modification 
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Figure  7.      Lookup  Table  After  Slope  Modification 
D.     APPLICATION  OF  LOOKUP  TABLE  OPERATIONS 
I.     Lookup  Table  Stretch  Routines 

As  shown  in  Chapter  II  of  this  thesis  (Figures  1  and  2),  the  histograms  for 
turbid  water  images  and  bottom  images  have  high  concentrations  of  pixels  in  the  center 
range  of  intensity  levels  and  few  or  zero  in  the  outer  ranges.  In  these  situations  contrast 
enhancement  may  be  achieved  by  using  lookup  table  stretch  methods  to  expand  the 
histogram  to  cover  the  full  range  of  available  intensity  levels.  [Ref.  5]  The  following 
experiment  was  performed.  The  histogram  (shown  in  Figure  8)  of  a  typical  image  from 
of  a  sequence  of  bottom  images  was  evaluated  visually  for  the  high  and  low  intensity 
values  of  the  high  pixel  concentration  area.  Then  a  linear  lookup  table  was  modified 
using  the  slope  varying  routine  described  in  this  chapter  so  that  the  top  intercept 
matched  the  high  end  of  intensity  level  concentration  and  the  bottom  intercept  matched 
the  low  end  of  intensity  level  concentration  as  shown  in  Figure  9.  The  histogram  routine 
was  then  applied  to  compare  the  enhanced  image  with  the  unenhanced  image.  The 
modified  histogram  is  shown  in  Figure  10.  The  result  was  an  image  whose  contrast  was 
significantly  enhanced.  The  problem  then  became  one  of  automating  the  histogram 
evaluation  process  and  making  it  responsive  to  changes  as  they  occurred  in  the  contin- 
uously acquired  images. 
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Figure  8.      Histogram  Before  Lookup  Table  Stretch 
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Figure  9.      The  Lookup  Table  Stretch 


15 


Figure  10.      Histogram  After  Lookup  Table  Stretch 

Histogram  equalization  provides  the  answer  to  the  histogram  evaluation  portion 
of  the  problem.   This  process  is  defined  by  a  transformation  function  of  the  form 


V(x) 


=       P(s)ds, 


where  p(s)  is  the  histogram  distribution  function,  5  is  a  dummy  variable  of  integration, 
x  is  a  normalized  input  gray  level,  and  V  is  the  transformed  output  gray  level.  This 
transformation  function,  when  plotted  for  the  full  gray  scale,  gives  the  shape  of  the  de- 
sired lookup  table.  [Ref  6]  The  histeq  function  uses  this  concept  to  find  the  output 
value  for  each  input  to  the  lookup  table  by  computing  the  normalized  cumulative  area 
up  to  that  input  intensity  level  on  the  histogram.  It  then  loads  those  output  values  into 
the  specified  lookup  table.  In  efTect,  the  function  steepens  the  slope  of  the  lookup  table 
plot  for  high  concentration  portions  of  the  histogram  and  shallows  the  slope  for  low 
concentration  portions,  a  localized  lookup  table  stretch.  Figures  11  and  12  show  the 
lookup  table  after  equalization  and  the  resulting  histogram.  Figures  13  and  14  show  an 
image  before  and  after  equalization.  Although  the  overall  image  appears  darker,  details 
of  the  bottom  are  more  readily  apparent. 
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Figure  11.      Lookup  Table  After  Equalization 


Figure   12.      Histogram  After  Equalization 
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Figure  13.      Image  Before  Equalization 


Figure   14.      Image  After  Equalization 
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The  problem  of  changing  image  frame  content  and  hence  changing  histogram 
was  solved  by  programming  two  alternatives  for  update  (performance  of  a  new 
histogram  equalization).  The  first  is  a  manual  update  performed  from  the  computer 
keyboard  by  the  operator  as  he  notices  the  image  content  changing  on  the  video  mon- 
itor. The  second  is  an  automatic  update  where  the  computer  routinely  grabs  a  frame 
and  performs  a  new  equalization  at  a  specified  time  interval.  With  either  alternative  the 
update  occurs  very  quickly  with  less  than  a  one  second  break  in  the  continuous  image 
acquisition  as  the  new  lookup  table  is  loaded. 

An  option  for  equalization  using  the  intensity  levels  of  a  specified  portion  of  the 
image  was  also  included.  The  user  designates  the  portion  of  the  image  using  the  same 
rectangle  described  in  the  histogram  portion  of  this  chapter.  The  rectangle  parameters 
are  used  to  define  the  area  used  for  the  equalization;  this  results  in  maximum  contrast 
for  that  area. 

2.     Other  Lookup  Table  Routines 

Chapter  II  of  this  thesis  also  describes  the  problems  caused  by  the  lighting  re- 
quired in  the  underwater  environment.  The  histogram  (shown  in  Figure  3)  of  an  image 
in  which  blooming  was  occurring  suggested  that  a  modification  to  the  portion  of  the 
lookup  table  containing  the  high  intensity  levels  was  in  order.  The  lookup  table  drawing 
routine  described  in  this  chapter  was  used  to  experiment  with  different  lookup  table 
shapes  which  reduced  the  intensity  of  high  value  pixels.  An  example  of  one  of  these  is 
shown  in  Figure  15.  Also,  the  invlut  and  abslut  functions  were  tried.  These  lookup  table 
shapes  are  shown  in  Figures  16  and  17.  The  procedures  most  effective  in  reducing  the 
blooming  without  affecting  the  remainder  of  the  image  used  the  inverse  and  absolute 
value  lookup  tables.   Both  were  included  as  menu  options. 

Finally,  the  loglut  function  was  tested  on  images  which  appeared  dark  due  to 
shadows.  The  logarithmic  values  had  a  very  definite  brightening  effect  on  the  images, 
similar  to  that  of  histogram  equalization  on  the  same  sequence.  The  loglut  function  was 
included  as  a  menu  option  as  a  simple,  one  keystroke  alternative  to  equalization  for  dark 
images.   The  logarithmic  lookup  table  is  shown  in  Figure  18. 

E.     APPLICATION  OF  OTHER  OPERATIONS 

Several  non-real  time  image  processing  functions  were  found  to  be  effective  in  en- 
hancing turbid  water  and  bottom  images  and  were  included  as  menu  options.  These 
functions  all  operate  on  a  single  image  frame  fixed  in  memory.  The  time  required  to 
perform  these  operations  varies  from  four  to  thirty  seconds  depending  on  the  operation 
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Figure  15. 
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Figure   16.      Inverse  Lookup  Table 
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Figure  17. 


Absolute  Value  Lookup  Table 
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Figure   18.      Logarithmic  Lookup  Table 
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and  computer  speed.  Analysis  of  the  effectiveness  of  these  options  was  limited  to  visual 
evaluation  of  the  actual  image. 

The  zoom  function  was  included  as  a  menu  option  to  assist  in  object  identification. 
In  order  to  allow  selection  of  the  portion  of  the  screen  to  be  enlarged,  a  rectangle  routine 
similar  to  that  used  in  the  histogram  and  equalization  options  was  programmed.  In  this 
case,  however,  the  rectangle  location  parameters  are  sent  to  the  pan  and  scroll  functions 
to  shift  the  image  in  memory  to  the  upper  left  quadrant,  the  portion  of  memory  used  in 
the  zoom  function.  Once  the  image  is  enlarged,  it  may  be  processed  by  other  options 
that  operate  on  a  single  frame  fixed  image. 

The  sharpen  function  was  included  as  an  option  to  provide  the  capability  to  amplify 
the  edges  of  a  fixed  image.  It  is  potentially  helpful  in  object  identification  or  in  finding 
an  entry  point  on  the  bottom. 

The  lopass  function  provides  the  capability  to  filter  out  high  frequecy  noise  from  a 
fixed  image.   It  is  included  as  a  menu  option  as  a  potential  tool  for  turbid  water  viewing. 

The  sobel  function  was  included  as  a  menu  option  to  provide  an  edge  detection  ca- 
pability.   It  is  potentially  useful  in  object  identification. 

Image  averaging  was  included  as  a  menu  option  using  the  average  function.  When 
this  option  is  called,  a  specified  number  of  image  frames  are  acquired  and  stored  to  disk. 
They  are  then  averaged,  pixel  by  pixel,  and  the  result  is  displayed  on  the  video  monitor. 
This  option  has  potential  usefulness  when  operating  in  turbid  water  with  the  camera 
mounted  on  a  stable  platform.  It  allows  the  selection  of  the  number  of  image  frames  to 
be  averaged  up  to  disk  memory  capacity.  It  is  the  slowest  of  the  processing  options  due 
to  the  high  number  of  read  write  to  disk  operations  and  large  number  of  calculations 
required. 
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V.     ON  SITE  TESTING 

A.  GENERAL 

An  important  aspect  of  this  project  was  to  provide  the  NUWES  recovery  section 
with  a  package  that  was  easy  to  use  and  helpful  in  their  specific  application.  For  this 
reason  a  dialogue  was  established  and  maintained  with  the  NUWES  engineers  and  re- 
covery section  operators  for  the  duration  of  the  development  process.  The  author  made 
three  visits  to  NUWES.  The  purpose  of  the  first  visit  was  to  observe  a  recovery  opera- 
tion first  hand  and  define  the  key  elements  of  the  problem.  The  next  visit  was  conducted 
midway  into  the  project.  Its  purpose  was  to  demonstrate  a  prototype  "Menu"  program 
and  obtain  feedback  about  the  usefulness  of  the  various  image  processing  options  and 
the  format  of  the  on-screen  presentation.  The  final  visit  was  an  on-site  test  of  the  pro- 
gram, conducted  on  a  torpedo  recovery  vessel  with  cameras  and  equipment  on  the  bot- 
tom under  actual  recovery  conditions.  The  aft  section  of  a  torpedo  was  used  as  the  test 
object  to  be  viewed.  The  results  of  that  test  and  the  feedback  obtained  from  NUWES 
operators  and  engineers  are  the  subject  of  this  chapter. 

B.  SPECIFIC  APPLICATIONS 

The  three  specific  applications  initially  identified  for  processing  were  turbid  water 
viewing,  bottom  viewing  while  searching,  and  blooming.  Additionally,  during  the  course 
of  the  development,  enhancement  of  shadowed  images  and  object  identification  became 
areas  of  interest. 

The  operations  expected  to  be  most  effective  during  turbid  water  viewing  were 
histogram  equalization,  the  low  pass  filter,  and  image  averaging.  This  portion  of  the  test 
provided  the  most  difficult  environment.  The  histogram  equalization  did  little  to  im- 
prove the  the  overall  quality  of  viewing.  The  periods  when  the  object  was  visible 
through  the  silt  were  so  short  that  by  the  time  a  manual  update  was  performed  the  ob- 
ject had  disappeared.  Similarly,  in  the  automatic  update  mode  the  randomness  of  the 
short  appearances  of  the  object  defeated  the  periodic  update  of  the  lookup  table.  The 
situation  was  further  complicated  by  a  horizontal  dark  gray  synchronization  bar  moving 
through  the  image  every  two  to  three  seconds.  This  reduced  the  level  of  effectiveness 
of  the  equalization  since  the  bar's  intensity  values  were  included  in  the  computation. 
In  fact,  the  equalization  process  enhanced  the  bar.  The  synchronization  bar  had  not 
been  visible  during  the  previous  on  board  visit  or  on  videotape  and  so  may  have  been  a 
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problem  peculiar  to  the  particular  equipment  used  on  this  test.  The  low  pass  filter  also 
proved  to  be  of  limited  value  in  this  environment.  Again,  the  periods  when  the  object 
was  visible  at  all  were  so  short  that  it  was  practically  impossible  to  freeze  an  image  to 
be  processed. 

Image  averaging  met  with  the  highest  degree  of  success  during  the  turbid  water 
portion  of  the  test.  With  the  recovery  device  (camera  platform)  sitting  firmly  on  the 
bottom,  the  image  average  routine  was  able  to  capture  those  fleeting  moments  when  the 
object  was  visible  along  with  the  time  when  it  was  not.  The  averaging  computation  then, 
to  some  extent,  removed  the  turbidity  from  the  image.  The  result  was  a  non-real  time 
image  in  which  the  object  was  recognizable.  The  time  elapsed  from  the  beginning  of  the 
process,  averaging  six  frames,  was  approximately  one  minute.  This  was  attributed  to  a 
relatively  low  computer  speed  (6  megahertz)  and  the  read  write  operations  to  the  hard 
disk.5 

The  operations  included  in  the  program  for  enhancement  of  images  during  bottom 
viewing  were  histogram  equalization  and  image  sharpening.  During  the  test  a  camera 
was  pointed  at  the  bottom  and  these  operations  were  performed.  The  histogram  equal- 
ization routine  produced  good  results  in  increasing  the  contrast,  making  bottom  features 
more  easily  identifiable.  The  synchronization  bar  again  reduced  effectiveness,  but  the 
value  of  the  equalization  was  proven  for  this  application.  The  sharpen  routine  was  also 
helpful  in  the  enhancement  of  bottom  images,  although  this  routine  does  not  operate  in 
real  time.  It  was  concluded  that  both  of  these  routines  would  be  of  considerable  value 
in  identifying  torpedo  entry  marks  on  the  bottom. 

The  inverse  and  absolute  value  lookup  table  operations  were  included  in  the  "Menu" 
program  to  help  solve  the  problem  of  blooming.  During  the  test  an  attempt  was  made 
to  induce  blooming  by  directing  the  lights  towards  the  test  object.  However,  due  to  the 
lighting  angle  and  relatively  dull  surface  of  the  object,  severe  blooming  as  seen  on  vide- 
otape did  not  occur.  The  inverse  and  absolute  value  operations  performed  as  expected 
on  the  minimal  blooming  induced,  and  operator  feedback  indicated  that  these  operations 
would  have  some  utility  during  an  actual  recovery.  Interest  was  also  expressed  in  the 
capability  to  experiment  with  different  lookup  table  shapes  using  the  lookup  table 
drawing  routine  and  save  those  which  helped  solve  the  blooming  problem. 


5  An  attempt  to  speed  up  the  image  averaging  operation  is  described  in  part  C  of  this  chapter. 
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The  zoom  operation  was  also  demonstrated  and  although  there  was  no  specific  ap- 
plication for  its  use  during  the  test,  it  was  agreed  that  it  would  be  helpful  in  certain  in- 
stances during  recovery  operations. 

C.     USER  FEEDBACK 

Upon  completion  of  the  test,  comments  were  solicited  from  the  recovery  section 
operators  and  engineers  concerning  ways  to  improve  the  "Menu"  program.  The  sug- 
gestion was  made  that  the  program's  menu  presentation  could  be  revised  to  allow  the 
use  of  a  mouse  to  make  selections  and  move  rectangles  and  pointers  on  the  video  mon- 
itor. A  high  degree  of  interest  was  also  expressed  in  improving  the  image  averaging 
routine.  Methods  suggested  included  compensating  for  motion  of  the  camera  platform, 
using  weighted  averaging  techniques,  and  speeding  up  the  process.  Once  ashore,  an  at- 
tempt to  provide  more  speed  was  made  by  installing  the  framegrabber  board  on  a  20 
megahertz  Compaq  386  computer  and  copying  the  "Menu"  program  to  a  virtual  disk. 
The  combination  of  higher  computer  speed  and  use  of  RAM  for  temporary  image  stor- 
age produced  a  significant  reduction  in  the  time  required  to  perform  the  operation  (ten 
frames  were  averaged  and  the  result  was  displayed  in  ten  seconds). 

Another  area  of  interest  was  the  potential  to  use  the  color  capability  available  by 
using  three  framegrabber  boards.  Suggestions  ranged  from  coloring  the  rectangle 
superimposed  on  image  frames  during  some  operations  to  creating  pseudocolor  images. 

In  general,  there  was  enthusiasm  expressed  for  the  use  and  continued  development 
of  the  capability  provided  by  this  image  processing  system. 
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VI.     CONCLUSIONS 

The  problem  of  enhancement  of  underwater  images  during  torpedo  recovery-  oper- 
ations was  analysed  and  broken  down  into  processing  for  three  specific  viewing  envi- 
ronments. These  were  viewing  through  turbid  water,  viewing  of  bottom  features,  and 
viewing  when  blooming  occurs.  The  tools  available  to  solve  these  problems  were  then 
studied.  These  included  an  IBM  AT  Personal  Computer  with  the  PCVISIONplus 
FRAMEGRABBER  installed  and  a  library  of  image  processing  functions  called  ITEX 
PCplus.  The  "Menu"  program  was  then  developed  by  incorporating  the  image  process- 
ing functions  into  an  interactive  program  with  a  menu  format.  The  program  provides 
convenient  and  easy  access  to  the  image  processing  functions  provided  by  ITEX  PCplus. 
It  also  provides  capability  for  automatic  or  manual  equalization,  a  segmentation  routine 
for  improved  contrast  on  a  selected  portion  of  an  image,  graphical  input  of  lookup  ta- 
bles, modification  of  existing  lookup  tables,  and  storage  of  lookup  tables  to  disk.  The 
histogram  and  lookup  table  display  routines  were  built  and  included  as  analysis  tools. 
These  provide  on  screen  display  of  lookup  tables  and  histograms  and  the  capability  to 
display  enhanced  and  unenhanced  versions  of  an  image  and  their  respective  histograms 
for  comparison  purposes.  Some  non-real  time  functions  included  were  image  averaging, 
a  Sobel  edge  detector,  a  low  pass  filter,  a  sharpening  routine,  and  a  zoom  feature  which 
allows  interactive  selection  of  the  portion  of  an  image  to  be  enlarged.  Documentation 
for  the  "Menu"  program  was  written  in  the  form  of  a  User's  Manual  (see  Appendix  A) 
which  includes  framegrabber  installation  instructions,  a  tutorial  on  the  use  of  the 
"Menu"  program,  and  a  quick  reference  list. 

A  test  of  the  final  package  was  conducted  aboard  a  recovery  vessel  at  XUWES  with 
cameras  and  equipment  on  the  bottom  under  actual  recovery  conditions.  The  program 
was  found  to  be  particularly  effective  when  viewing  bottom  features  and  during  bloom- 
ing. The  equalization  procedures,  lookup  table  drawing  capability,  absolute  value 
lookup  table  modification,  image  averaging,  and  the  zoom  function  were  all  determined 
to  be  potentially  valuable  tools  for  actual  recovery  operations.  Appreciation  was  also 
expressed  for  the  documentation  provided  with  the  program. 

More  development  was  called  for  in  the  turbid  water  viewing  environment,  specif- 
ically with  regard  to  the  improvement  of  the  image  averaging  routine.  The  current 
function  averages  a  finite  number  of  image  frames  chosen  by  the  user,  each  weighted 
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equally.  All  the  frames  are  stored  on  the  hard  disk  and  then  read  back  into  RAM  for 
averaging.  The  speed  of  the  process  is  limited  by  the  clock  speed  of  the  computer  and 
the  read  writes  to  and  from  disk.  A  stable  camera  platform  is  essential  or  the  resulting 
image  becomes  extremely  blurred.  Improvements  could  include  continuous  weighted 
frame  averaging,  compensating  for  camera  motion  by  cross-correlating  image  frames, 
and  use  of  a  faster  computer  with  a  virtual  disk  for  increased  speed. 

The  "Menu"  program  with  documentation  was  delivered  to  NUWES  for  operational 
use. 
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APPENDIX  A.     USER'S  MANUAL 

A.  GENERAL 

The  "Menu"  software  package  is  a  menu-driven  program  which  provides  real  time 
continuous  video  image  processing  using  lookup  table  modification  techniques.  It  also 
provides  limited  non-real  time  single  frame  enhancement  capability.  Control  of  the 
routines  in  "Menu"  varies  from  manual  (requiring  operator  input  to  make  an  adjust- 
ment) to  automatic  or  "hands-ofT  operation.  The  processing  routines  were  chosen  to 
provide  maximum  effectiveness  in  situations  where  the  images'  pixel  values  are  limited 
to  a  relatively  small  region  of  the  total  gray  scale  range.  Typically  this  would  include 
images  of  objects  which  have  actual  small  variations  in  shading  or  those  whose  contrast 
has  been  reduced  by  obscuration  or  dim  lighting. 

B.  EQUIPMENT  REQUIRED 

The  following  equipment  is  required  to  use  the  "Menu"  program: 

•  IBM  Personal  Computer  AT,  XT,  or  PC6  or  100%  compatibles  with  Extended 
Graphics  Adaptor  card  and  monitor. 

•  Imaging  Technology  PCVISIONplus  Frame  Grabber"?  card  with  standard  cable. 
(NOTE:   The  older  model  PCVISIOX  Frame  Grabber  will  not  work.) 

•  Analog  video  monitor  (such  as  SONY  Trinitrons  Model  1271Q). 

•  Video  source  (camera  or  video  tape  player). 

C.  SIMPLE  SETUP 

The  PCVISIONplus  Frame  Grabber  User's  Manual  provides  all  the  information 
required  to  place  the  Frame  Grabber  into  operation.   Here  is  a  much  condensed  version 
of  that  procedure  with  some  helpful  hints  to  ease  the  way. 
1.     Getting  the  Board  Ready 

Figure  2-1  on  page  2-2  of  the  Frame  Grabber  Manual  shows  the  jumper  lo- 
cations on  the  board  and  Table  2-1  on  page  2-3  tells  how  they  are  configured  from  the 
factor}'.   There  are  three  position  functions  to  be  concerned  with: 
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•  Since  an  EGA  card  is  needed  to  run  "Menu"  the  memory  base  address  must  be 
changed  on  the  Frame  Grabber  board.  Address  D0000  is  recommended  since  this 
is  the  default  value  in  the  software.  To  change  to  address  D0000  simply  insert  a 
jumper  at  location  J8.  If  other  than  D0000  is  chosen  it  should  be  noted  for  entry 
each  time  "Menu"  is  run. 

•  "Menu"  has  a  default  register  base  address  of  100.  To  change  the  FrameGrabber 
board  insert  a  jumper  at  location  Jll.  Again,  if  100  is  not  used,  the  chosen  value 
will  have  to  be  entered  each  time  "Menu"  is  run. 

•  If  an  AT  host  is  being  used  a  jumper  is  required  at  location  J20.  If  an  XT  or  PC 
is  being  used  no  jumper  should  be  installed  at  J20. 

2.  Installing  the  Board 

Pages  3-1  through  3-3  of  the  FrameGrabber  Manual  provide  a  detailed  de- 
scription for  installing  the  board.  A  summary  is  as  follows: 

•  Turn  off  and  unplug  the  computer  and  periphials. 

•  Remove  the  cover. 

•  Place  the  board  in  any  open  slot  where  it  fits  and  install  screws. 

•  Reinstall  cover. 

•  Plug  in  and  turn  on  computer. 

3.  Cable  Installation 

Page  3-5  of  the  Frame  Grabber  Manual  shows  how  to  connect  the  cable.  Here's 
the  bottom  line: 

•  Plug  the  rectangular  female  connector  on  the  cable  into  the  top  connector 
protruding  from  the  board. 

•  Connect  the  white  labeled  BXC  connector  to  the  video  source  output. 

•  Connect  the  green  labeled  BNC  connector  to  the  'line'  input  (not  the  'green'  input 
of  the  red,  green,  blue  connections). 

•  Do  not  connect  the  red-  or  blue-labeled  cables. 

4.  Running  the  Diagnostics. 

Page  3-10  of  the  Frame  Grabber  Manual  describes  the  procedure  for  running 
the  diagnostic  software  included  with  the  Frame  Grabber.  Since  the  factory  configura- 
tion has  been  changed  as  described  above,  the  command  <  configure  >  will  have  to  be 
entered  followed  by  the  register  base  address  (100)  and  the  memory  base  address 
(D0000).  The  diagnostics  take  several  minutes  to  run,  so  be  patient.  As  stated  in  the 
manual,  system  memory  or  other  devices  mapped  into  the  same  memory  space  or  regis- 
ter addresses  will  cause  some  tests  to  fail. 
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D.     LET'S  TRY  IT  OUT 

This  is  an  introduction  to  the  "Menu"  program.  It  lets  you  walk  through  some  of 
the  features  of  the  program;  sort  of  a  hands-on  familiarization.  Turn  on  your  source 
(camera  or  VCR)  and  video  monitor  now  to  let  them  warm  up.  On  your  computer  go 
to  the  directory  where  you  want  to  keep  the  program  and  copy  'menu.exe'  there. 

1.  Load  the  Program 

To  get  started  simply  type  the  word  'menu'  followed  by  <  enter  >.  No  <  en- 
ter >  is  needed  for  any  selection  from  here  on  out  unless  otherwise  noted.  The  menu 
should  now  be  displayed  on  the  computer  monitor.  If  you  are  not  using  the  memory 
base  address  (D0000)  and  register  base  address  (100)  as  described  in  paragraph  C  above, 
hit  'K'  to  change  setup  and  follow  the  instructions  on  your  monitor. 

2.  Live  Action 

Menu  selection  E'  provides  continuous  acquisition  and  display  of  the  incoming 
real  time  video  images.  Hit  'E'  and  you  should  see  live  action  from  your  camera  or  VCR 
(if  you  are  using  a  tape,  make  sure  you  press  the  play  button  first). 

3.  Take  a  Snapshot 

Now  let's  freeze  the  action.  Hit  'F'  for  snap.  This  acquires  and  displays  a  single 
frame.  You  should  see  a  stop  action  snapshot  on  your  video  monitor.  The  image  is  also 
stored  now  in  framegrabber  memory  A,  the  default  memory  chosen  when  you  started 
up. 

4.  Take  Two  Snapshots 

Now  hit  'D'  for  display  memory  B.  Take  another  snapshot  with  the  'F'  key. 
The  second  snapshot  is  now  displayed  on  the  monitor  and  stored  in  framegrabber 
memory  B.  Try  going  back  to  memory  A  (the  C  key)  to  make  sure  the  first  snapshot 
is  still  there.  Note  that  we  are  using  the  display  memory  keys  and  not  the  select  memory 
keys  to  perform  these  operations.  More  on  select  memory  later  when  we  start  storing 
images. 

5.  All  Black 

Let's  erase  one  of  our  snapshots.  Display  memory  B  (the  D'  key)  and  then  hit 
the  'G'  key  for  screenclear(O).  This  clears  the  screen  and  framegrabber  memory  B  to 
black  or  intensity  level  zero. 

6.  All  White 

While  still  in  memory  B  hit  the  'H'  key  for  screenclear(255).  This  erases  the 
screen  and  memory  B  to  white  or  intensity  level  255.  Darker  gray  shades  have  lower 
intensity  levels  while  lighter  shades  have  higher  intensity  levels. 
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7.  Put  Your  Snapshot  in  the  Photo  Album 

Saving  an  image  to  disk  is  analogous  to  putting  your  snapshot  in  a  photo  al- 
bum. You  should  still  have  an  image  in  frame  memory  A.  If  not,  use  the  display  mem 
A  and  snap  routines  to  get  one.  Now  hit  the  'A'  key  to  select  memory  A.  This  selects 
memory  A  for  data  manipulations  or  storage.  For  our  purposes  we  have  selected  the 
image  data  in  memory  A  to  store  on  disk.  Now  hit  the  T  key  and  you  will  be  asked  to 
enter  an  image  name.  Type  a  name  (e.g.,  'image')  to  store  the  image  on  the  current  disk 
and  directory  or  path  and  name  (e.g.,  'a:image')  to  store  it  in  a  different  place.  Then 
you'll  be  prompted  for  a  comment.  You  may  give  it  a  label  or  simply  type  'none'  (note 
that  <  enter  >  is  required  for  each  of  these  inputs).  The  disk  drive  selected  should  make 
some  noise  and  your  snapshot  is  safely  tucked  away  in  the  "photo  album"  (i.e.,  the  image 
file  is  written  to  the  disk). 

8.  Open  the  Album  to  See  Your  Snapshot 

Let's  see  if  the  image  was  really  saved.  First  erase  both  frame  memories  as 
described  in  paragraph  5  or  6  above.  Then  select  a  memory  using  'A'  or  B'.  Now  hit 
the  'J'  key.  Type  the  image  file  name  (followed  by  <enter>)  to  read  the  image  from 
disk.  The  image  should  appear  on  the  video  monitor.  If  it  does  not,  display  the  memory 
you  selected  (the  C  or  'D'  key)  and  the  image  will  appear. 

9.  Simple  Processing 

Now  let's  process  real  time  video  images  using  some  simple  lookup  table  mod- 
ification techniques.  Hit  the  'V  key  -  its  called  analysis,  but  really  it  just  lets  you  see 
how  your  image  is  being  changed.  The  graph  you  see  on  your  computer  monitor  is  the 
current  lookup  table.  The  horizontal  coordinate  represents  input  pixel  values  and  the 
vertical  represents  output  values. 

On  the  top  of  the  computer  monitor  screen  are  the  current  options  available  for 
selection.  The  snap  ('F')  and  grab  ('G')  options  work  the  same  as  on  the  previous  menu. 
Try  them!  Now  try  the  invlut  (T)  option  and  see  how  the  lookup  table  changes  (the 
blue  line  is  the  current  lookup  table  and  the  green  is  a  reference  linear  lookup  table  line). 
The  effect  on  the  image  is  just  like  looking  at  the  negatives  of  a  photograph,  except  that 
we  still  have  live  action.  Now  hit  the  L'  key  for  linear  lookup  table.  We  are  back  to 
unprocessed  live  images. 

Keys  'P'  (logarithmic)  and  'A'  (absolute  value)  also  modify  the  lookup  table. 
You  can  experiment  with  these  functions  to  see  the  effect  on  the  real  time  image  and 
their  corresponding  lookup  table  shape. 
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10.  Histogram  Equalization 

Now  for  the  most  powerful  of  our  lookup  table  modification  routines.  Hit  'E' 
for  equalize.  You  now  are  given  two  choices,  default  and  rectangle.  This  determines 
which  portion  of  the  image  will  be  used  to  do  the  processing.  Most  of  the  time  you 
will  probably  use  the  default  mode  but  let's  try  rectangle  now  just  for  fun.  Hit  'R'  for 
rectangle.  On  the  video  monitor  you  will  see  the  live  images  being  displayed  with  a 
rectangle  overlaid.  Now  hit  'Num  Lock'  on  your  keyboard.  Using  the  '2',  '4',  '6',  and 
'8'  keys  on  the  number  pad  (the  ones  with  the  arrows)  you  can  move  the  rectangle 
around  on  the  screen.  The  corner  numbers  (T,  '3',  '7',  and'9')  are  used  to  change  the 
size  of  the  rectangle.  If  your  image  had  a  large  white  or  black  area  in  which  you  were 
not  interested,  it  would  lessen  the  impact  of  the  equalization  routine  on  the  area  of  in- 
terest. So  you  would  frame  the  area  of  interest  with  the  rectangle.  After  you  have 
placed  the  rectangle  where  you  want  it,  hit  S'  for  stop. 

Now  you  have  to  choose  between  manual  or  automatic  update.  Try  manual 
update  first.  Hit  the  'M'  key.  The  image  will  equalize,  the  new  lookup  table  is  plotted, 
and  live  action  continues  using  the  new  lookup  table.  In  this  mode,  as  you  notice  the 
contents  of  the  image  changing  you  can  hit  the  'U'  key  and  equalization  will  occur  again. 

Let's  go  back  now  and  try  the  default  rectangle  parameters.  Hit  C  to  change 
parameters.  Now  hit  'D'  for  default  and  'C  for  continuous  update.  You  are  prompted 
for  time  between  updates  in  seconds.  Enter  any  number  2  through  9.  The  computer 
will  perform  a  new  equalization  at  the  interval  you  specify.  This  is  intended  for  situ- 
ations where  the  contents  of  the  image  are  changing  often  and  you  don't  want  to  have 
to  remember  to  update.   When  you  have  seen  enough  updates  hit  'S'  for  stop. 

11.  Before  and  After 

Let's  look  at  some  histograms  and  see  whats  happening  when  we  equalize.  Hit 
the  'H'  key  for  histogram.  Again,  you  may  choose  between  default  or  rectangle  to 
specify  a  portion  of  the  image.  Hit  'D'  for  default.  Now  hit  the  'L'  key  for  computation 
of  linear  lookup  table  histogram.  After  the  "COMPUTING  HISTOGRAM"  indicator 
goes  away,  hit  'M'  for  computation  of  modified  histogram.  When  the  indicator  goes 
away  this  time,  hit  'D'  to  display  the  linear  lookup  table  histogram.  This  is  the  infor- 
mation used  to  perform  the  equalization.  The  image  on  the  video  monitor  is  the  un- 
processed version.  Now  hit  'E'  to  display  the  modified  histogram.  The  video  image  is 
enhanced  and  the  "spread  out"  histogram  is  displayed.  As  you  toggle  back  and  forth 
between  'D'  and  'E'.  you  can  see  graphically  how  equalization  lets  us  use  the  full  range 
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of  intensity  levels  and  the  corresponding  enhancement  of  contrast  on  the  image.    Con- 
vinced?  Hit  the  'S'  key  and  let's  try  something  else. 

12.  Modifying  the  Lookup  Table 

Hit  the  'N'  key.  Now  we'll  go  again  to  the  number  pad  with  'Num  Lock'  on. 
Using  the  corner  keys  (T,  '3',  '7',  and  '9')  we'll  change  the  slope  of  the  lookup  table  by 
moving  the  top  or  bottom  of  the  curve.  To  shallow  out  the  lookup  table  use  the  '9'  key- 
to  move  the  top  to  the  right  and  the  T  key  to  move  the  bottom  to  the  left.  Notice  that 
enhancement  to  the  live  images  changes  as  you  modify  the  curve.  When  you  are  finished 
hit  'S'  to  stop  and  the  resulting  lookup  table  is  shown.  If  you  were  to  develop  a  "fa- 
vorite" lookup  table  which  you  think  you  might  want  to  use  again,  there  is  a  way  to  save 
it  in  the  computer.   We'll  talk  about  that  later. 

13.  Draw  Your  Own  Lookup  Table 

Now  hit  'M'  to  perform  another  form  of  manual  lookup  table  modification. 
We'll  do  this  by  drawing  lines.  Using  the  number  pad  with  'Num  Lock'  on  move  the 
pointer  on  the  computer  monitor  by  hitting  the  '2',  '4',  '6',  and  '8'  keys.  This  will  move 
it  one  pixel  in  the  direction  of  the  arrows  on  the  keys.  The  next  key  clockwise  from  each 
of  these  will  move  it  in  the  same  direction  ten  pixels  .  (E.g.,  The  T  key  moves  the 
pointer  in  the  same  direction  as  the  '2'  key  only  faster.)  When  the  pointer  is  at  a  lo- 
cation where  you  want  the  line  to  begin  hit  'B'  for  begin  line.  Then  move  the  pointer 
to  a  location  where  you  want  the  line  to  end  and  hit  'E'  for  end  line.  Now  hit  'D'  to 
draw  the  line  and  your  change  will  take  effect  on  your  image.  You  may  draw  other  lines 
to  connect  segments  and  build  a  whole  new  lookup  table.  To  get  rid  of  the  clutter  just 
hit  S'  and  then  reenter  the  routine  by  hitting  'M'. 

Now  hit  'P'  to  put  a  pointer  on  the  live  video  image.  To  move  the  pointer  use 
the  number  pad  keys  in  the  same  manner  as  you  did  to  move  the  pointer  on  the  com- 
puter monitor  described  above.  Notice  on  your  computer  monitor  that  the  position  of 
the  pointer  and  the  intensity  level  are  displayed.  One  possible  use  of  this  capability 
would  be  to  find  the  intensity  level  of  an  object  in  the  image  and  use  the  draw  lookup 
table  routine  to  change  all  pixels  with  that  intensity  to  another  level.  Again,  if  you  were 
to  draw  a  really  good  all-purpose  lookup  table,  you  could  save  it  with  a  routine  which 
we'll  describe  later.  Hit  'S'  to  exit  the  video  monitor  pointer  routine  and  'S'  again  to  exit 
the  modifv  routine. 
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Finally,  let's  go  back  to  the  main  menu  to  see  the  lookup  table  storage/selection 
routines  and  try  some  single  frame  (frozen)  processing.  Type  'S'  to  exit  the  analysis 
function. 

14.  Storing  and  Retrieving  Your  Favorite  Lookup  Tables 

Remember  that  lookup  table  you  designed  using  the  modify  routine  a  while  ago? 
Here's  how  you  can  save  it  for  future  use.  Hit  the  'M'  key.  Type  in  a  name  (e.g.  'mylut' 
or  'a:mylut'),  hit  enter,  and  the  values  are  saved.  To  get  it  back  hit  the  'N'  key  to  read 
the  lookup  table  values  from  disk.  Again,  type  in  the  File  name  followed  by  enter  and 
your  lookup  table  will  take  effect. 

15.  Piggyback  Lookup  Tables 

The  framegrabber  board  has  two  sets  of  lookup  tables,  one  as  the  image  is  put 
into  frame  memory  and  one  as  the  image  is  output  to  the  monitor.  Up  to  this  point  we 
have  only  used  the  output  lookup  table.  Now  we'll  modify  both  sets  to  achieve  a  dif- 
ferent result.  Hit  the  'P'  key  to  select  the  input  lookup  table.  Then  go  back  to  the 
analysis  routine  ('L')  and  do  a  histogram  equalization.  Everything  still  looks  the  same 
as  the  last  time  we  equalized,  including  the  video  image,  but  this  time  the  computed 
values  went  into  the  input  lookup  table.  Now  exit  the  equalization  and  analysis  routines 
to  get  back  to  the  main  menu.  Hit  the  'O'  key  to  select  the  output  lookup  table.  Then 
go  back  into  analysis.  You  may  now  use  any  of  the  lookup  table  modification  routines 
(e.g.  invlut,  abslut,  etc.)  and  you  will  get  a  "piggybacked"  lookup  table  effect  on  the  im- 
age. To  get  back  to  normal,  hit  L'  to  linearize  the  output  lookup  table.  Go  back  to  the 
main  menu  and  select  and  linearize  the  input  lookup  table.  Note  that  when  you  use  this 
procedure,  always  put  values  derived  by  equalization  into  the  input  lookup  table  or  the 
equalization  will  be  based  on  already  modified  data. 

16.  Zoom 

Now  hit  the  'R'  key  to  start  the  zoom  function.  On  the  video  monitor  you  can 
see  the  rectangle  we  saw  before.  Using  the  number  pad  with  "Sum  Lock'  on  you  can 
move  the  rectangle  with  the  '2',  '4',  '6',  and  '8'  keys  in  the  direction  of  the  arrows  on  the 
keys.  The  size  of  this  one  is  fixed  at  approximately  one  quarter  of  the  screen,  since  we 
will  be  expanding  the  image  to  four  times  the  original.  Move  the  rectangle  to  the  por- 
tion of  the  image  you  wish  to  zoom  in  on  and  hit  'S'.  The  expanded  image  is  displayed. 
Notice  that  some  resolution  is  lost  from  the  original  because  we  are  "zooming"  digitally 
as  opposed  to  changing  the  focal  length  of  the  camera  lens.  Hit  'M'  to  get  back  to  the 
menu.  Now  you  can  process  the  "zoomed"  image  using  the  fixed  image  processing  rou- 
tines described  below. 
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17.  Unzoom 

To  return  to  a  normal  image  after  using  the  zoom  routine  and  performing  any 
fixed  image  processing,  hit  the  'R'  key.  It  is  essential  that  this  be  done  when  you  are 
finished  with  a  "zoomed"  image  or  the  zoom  parameters  will  interfere  with  other  rou- 
tines. 

18.  Sharpen 

Snap  an  image  by  hitting  the  'F'  key  or  by  using  the  zoom  routine  (the  'R'  key). 
Then  hit  'S'  to  sharpen  it.  If  you  watch  closely  you  can  see  the  sharpening  occurring 
on  the  displayed  image. 

19.  Low  Pass  Filter 

Snap  or  zoom  another  image  and  then  hit  "T  to  call  the  low  pass  filter  routine. 
This  function  is  designed  to  remove  high  frequency  noise  from  a  fixed  image.  Again, 
looking  closely  you  can  see  the  change  occurring. 

20.  An  Edge  Detector 

One  more  time,  snap  or  zoom  an  image.  This  time  hit  the  X"  key  to  perform 
a  Sobel  Edge  Filter.  Note  the  interesting  display  of  the  image's  edges  on  the  video 
monitor. 

21.  Image  Averaging 

Hit  the  'V  key.  Then  enter  the  number  of  image  frames  you  want  to  average 
together.  This  number  may  be  limited  by  the  amount  of  disk  space  you  have  available 
(2  is  a  good  number  to  start  with).  When  the  average  image  is  displayed,  you  should 
notice  that  objects  which  passed  across  the  field  of  view  during  the  operation  have  faded 
away.  The  more  frames  you  average,  the  less  visible  these  objects  will  become.  Image 
averaging  is  very  effective  at  eliminating  transient  noise,  but  requires  a  stable  platform 
for  the  source. 

To  exit  the  program  hit  the  'W  key. 

E.     QUICK  REFERENCE  LIST 

This  list  presents  the  program's  options  in  the  order  they  appear  with  a  brief  de- 
scription of  their  output. 
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1.  A  for  select  mem  A 

Selects  framegrabber  memory  A  for  storing  or  reading  an  image  to/from  disk. 

2.  B  for  select  mem  B 

Selects  framegrabber  memory  B  for  storing  or  reading  an  image  to/from  disk. 

3.  C  for  display  mem  A 

Selects  framegrabber  memory  A  for  acquisition  of  image  from  camera  and  dis- 
plays memory  A  on  the  monitor. 

4.  D  for  display  mem  B 

Selects  framegrabber  memory  B  for  acquisition  of  image  from  camera  and  dis- 
plays memory  B  on  the  monitor. 

5.  E  for  grab 

Continuously  acquires  image  frames  from  source  and  displays  on  video  monitor. 
Uses  framegrabber  memory  designated  by  'display  mem'  function. 

6.  F  for  snap 

Acquires  a  single  frame  from  source  and  displays  it  on  the  video  monitor.   Uses 
framegrabber  memory'  designated  by  'display  mem'  function. 

7.  G  for  sclear(O) 

Clears  screen  and  memory  to  intensity  level  0. 

8.  H  for  sclear(255) 

Clears  screen  and  memory  to  intensity  level  255. 

9.  I  to  store  image 

Stores  contents  of  selected  memory  to  disk.    Prompts  input  of  filename  (includ- 
ing path  if  desired)  and  comment. 

10.  J  to  read  image 

Reads  image  from  disk  file  to  selected  framegrabber  memory.    Prompts  input 
of  filename  (including  path  if  desired). 

11.  K  to  change  setup 

Allows  change  of  software  memory  base  address  and  register  base  address.    On 
screen  instructions  and  example  inputs  provided. 

12.  L  for  analysis 

Displays  lookup  tables  and  histograms.  Performs  equalization  and  lookup  table 
modification  routines.   Options  as  follows: 
a.     m  =  modlut 

Modifies  lookup  table  by  drawing  lines.    Places  pointer  on  video  monitor. 
Move  pointer  to  desired  start  point  and  end  point  using  'Num  Lock'  on  and  number  pad 
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keys  '2',  '4',  '6',  and  '8'  for  down,  left,  up,  and  right  respectively.    Use  keys  T,  '!',  '9', 
and'3'  for  faster  movement  down,  left,  up,  and  right  respectively. 

(1)  b=beginline.     Marks  beginning  of  line. 

(2)  e=endline.     Marks  end  of  line. 

(3)  d=  drawline.     Draws  line  and  enters  new  values  into  lookup  table." 

(4)  p  =  pointer.  Places  pointer  on  video  monitor.  Gives  pointer  location 
and  pixel  intensity  level  of  that  location  on  computer  monitor.  Use  same  method  to 
move  pointer  as  for  line  drawing  routine  above. 

(5)  5=  stop.     Exit  lookup  table  modify  routine. 

b.  n  =  mod2  lut 

Changes  slope  of  existing  lookup  table  and  implements  results.  With  'Num 
Lock'  on,  use  number  pad  keys  T  and  '3'  to  move  bottom  of  curve  left  and  right  re- 
spectively. Use  '7'  and  '9'  keys  to  move  top  of  curve  left  and  right  respectively.  Hit  s' 
to  exit  this  routine. 

c.  h  =  histogram 

Computes  and  displays  histograms  for  images  processed  with  linear  and 
modified  lookup  tables.  Snaps  images  before  and  after  processing  and  allows  for  com- 
parison. 

H  j  d=  default.  Selects  default  portion  of  image  for  histogram  computa- 
tion (approximately  90%  of  image). 

(2)  r=  rectangle.  Displays  rectangle  superimposed  on  image  on  video 
monitor.  Change  portion  of  image  used  for  histogram  computation  by  altering  rectangle 
size  and  location.  With  'Num  Lock'  on,  use  number  pad  keys  '2\  '4',  '6'.  and  '8'  to  move 
rectangle  down,  left,  right,  and  up  respectively.  Use  '1'  and  '3'  keys  to  increase  and 
decrease  height.  Use  '7'  and  '9'  keys  to  increase  and  decrease  width.  Hit  's'  to  exit  rec- 
tangle routine. 

f3j  1=  com  linhist.  Computes  histogram  of  image  processed  by  linear 
LUT.   Displays  "COMPUTING  HISTOGRAM"  message  until  done. 

(4)  m=  com  modhist.  Computes  histogram  of  image  processed  by 
modified  LUT.   Displays  "COMPUTING  HISTOGRAM"  message   until  done. 

(5 j  d=dis  linhist.  Displays  histogram  computed  above  using  image 
processed  by  linear  lookup  table.   Unprocessed  image  shown  on  video  monitor. 

(6 j  e  =  dis  modhist.  Displays  histogram  computed  above  using  image 
processed  by  modified  lookup  table.   Processed  image  shown  on  video  monitor. 

f7 j     s=  stop.     Exit  histogram  routine. 
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d.  e=  equalize 

Performs  histogram  equalization  with  update  on  command  or  automatically 
at  specified  intervals.   Displays  resulting  lookup  table  after  each  update. 

(1)  d=  default.  Selects  default  portion  of  image  to  use  for  equalization 
(approximately  90%  of  image). 

(2)  r=  rectangle.  Allows  use  of  rectangle  to  specify  portion  of  image  to 
be  used  as  described  in  histogram  rectangle  routine  above  (paragraph  12.c.(2)). 

(3)  c  =  continuous  update.  Selects  automatic  update  of  lookup  table  at 
specified  intervals.  Prompts  for  desired  interval.  Hit 's'  to  exit  equalization  routine  from 
this  mode. 

(4)  m  =  manual  update.     Selects  lookup  table  update  on  command. 

(a)  u  =  update—    Performs  new  equalization  and  updates  lookup 
table. 

(b)  c  =  changeparam—     Returns    to    default/rectangle   portion    of 
image  options  for  changing  parameters. 

(c)  f=  freezeframe—    Acquires  and  displays  a  single  frame  (snap). 

(d)  g=grab—    Continuous    acquisition    and    display    of   image 
frames. 

(e)  s=  stop—     Exit  equalization  routine  from  this  mode. 

e.  l=linlut 

Sets  lookup  table  to  linear  values.   Results  in  unprocessed  image. 

f.  i=invlut 

Sets  lookup  table  to  inverse  values.   Produces  a  negative  image. 

g.  p  =  loglut 

Sets  lookup  table  to  logarithmic  values.    Increases  image  brightness. 
h.     a=  abslut 

Inverses  intensity  levels  greater  than  127.   Reduces  image  brightness. 
i.    g  =  grab 

Continuous  aquisition  and  display  of  image  frames. 
j.    f=  freeze 

Aquisition  and  display  of  a  single  frame  (snap). 
A.     s  = stop 

Exit  analvsis  routine. 
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13.  M  to  store  lookup  table 

Stores  contents  of  selected  lookup  table  to  disk.  Prompts  input  of  filename 
(including  path  if  desired). 

14.  N  to  read  lookup  table 

Reads  lookup  table  values  from  disk  to  lookup  table  they  were  stored  from. 
Prompts  input  of  filename  (including  path  if  desired). 

15.  O  to  select  output  lookup  table 

Selects  the  output  lookup  table  to  perform  modifications  on  and  for  display  in 
the  analysis  routine.  When  used  in  conjunction  with  'P'  (select  input  lookup  table)  al- 
lows for  processing  through  two  modified  lookup  tables. 

16.  P  to  select  input  lookup  table 

Selects  the  input  lookup  table  to  perform  modifications  on  and  for  display  in  the 
analysis  routine.  When  used  in  conjunction  with  'O'  (select  output  lookup  table)  allows 
for  processing  through  two  modified  lookup  tables. 

17.  Q  for  zoom 

Enlarges  chosen  area  of  image.  Provides  moveable  rectangle  to  select  area  for 
zoom.  With  "Xum  Lock"  on  use  number  pad  keys  '2',  'A',  '6',  and  '8'  to  move  rectangle 
down,  left,  right,  and  up  respectively.  Hit  z'  to  stop  rectangle  and  fix/zoom  image.  Hit 
'm'  to  return  to  main  menu. 

18.  R  for  unzoom 

Returns  image  to  normal  state. 

19.  S  for  sharpen 

Sharpens  fixed  image  using  a  high  pass  filter.    May  be  used  on  zoomed  image. 

20.  T  for  lo»  pass 

Performs  a  low  pass  filter  operation  on  a  fixed  image.  Reduces  high  frequency 
noise.  Blurs  the  image.    May  be  used  on  zoomed  image. 

21.  U  for  edge  detector 

Performs  Sobel  edge  detector  operation  on  a  fixed  image.  Displays  detected 
edges.   May  be  used  on  zoomed  image. 

22.  V  for  image  average 

Averages  specified  number  of  image  frames  and  displays  results.  Prompts  input 
of  number  of  frames  to  average. 

23.  W  for  framegrabber  off  and  exit  program 
Terminates  program. 
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APPENDIX  B.     PROGRAM  LISTING 

/*This  program  provides  image  processing  options  in  a  menu  format.   It 
was  written  to  be  compiled  using  Microsoft  C  5.0.   It  also  uses  include 
files  and  the  library  of  functions  from  the  ITEX  PCplus  Image  Processing 
Package.   The  Imaging  Technology  PCplus  Framegrabber  must  be  installed 
to  run  the  program.   Graphics  in  the  program  require  an  EGA  card  and 
monitor. */ 

/*List  include  files*/ 

//include  <itexpfg.  h>  /*from  ITEX  PC  plus*/ 

//include  <stdtyp.  h>   /*from  ITEX  PCplus*/ 

//include  <stdio.  h> 

//include  <stdlib.  h> 

//include  <process.h> 

//include  <time.  h> 

//include  <conio.  h> 

//include  <graph.  h> 

/*Declare  variables*/ 

int 

i,err,n,m,t[ 10] ,s[ 10] , j ,highcut ,r, array [ 256] ,pp,qq,intens,c,a,b,da,db, 

e, co lor ,xl,x2,x3,x4,yyl,y2,y3,y4,mmm,nnn,res, 1, low, high, blx,bly,elx,ely, 

ppl,qql,d, lowcut , start , range, int ensity,y,g,h, array 2  [ 256] ,vt ,num; 

unsigned  regbase; 

unsigned  int  grp; 

unsigned  long  membase; 

char  name[  20]  ,  comment [  200]  ,  lutname[  20]  ,  lutname2[  20]  ; 

short  grmode; 

long  histvalsl[  256] ;  histvals2[  256]  ;  histvals[  256]  ; 

FILE  ''-stream; 

t  ime_t  1 1  ime , xx; 

float  dely ,delx, f,af,bf,gf,hf ,yi, delta; 

struct  videoconfig  config; 

/*Main  program*/ 
main( ) 


/*Set  up  the  framegrabber  board*/ 
sethdw( 0x100 , 0xD0000L,DUAL) ; 
setdim(512,512,8); 
fgon(); 
initialize( ); 
setvmask(OxOO); 
sclear(255); 

/*Initial  rectangle  and  LUT  parameters*/ 
a=30; b=30; da=450; db=450; grp=GREEN; 
loop: 

/*Present  menu  on  screen*/ 

_clearscreen(_GCLEARSCREEN); 
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printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 
printf 


"type:\n"); 

"A  for  select  mem  A\n"); 

MB  for  select  mem  B\n"^); 

"C  for  display  mem  A\n'); 

"D  for  display  mem  B\n"); 

"E  for  grab\n,r); 

"F  for  snap\n"); 

"G  for  sclear(0)\n"): 

"H  for  sclear(255)\n); 

"I  to  store  image\n"); 

11 J  to  read  image\n"): 

"K  to  change  setup\n'); 

"L  analysis\nM); 

"M  for  store  lut\nM); 

"N  for  read  lut\n"); 

"0  for  select  output  lut\n"); 

"P  for  select  input  lut\nM); 

"Q  for  zoom\n"): 

"R  for  unzoom\n'<); 

"S  for  sharpen\n '"); 

"T  for  low  pass\n'); 

"U  for  edge  detector\nM); 

"V  for  image  average\n"); 

"W  for  framegrabber  off  and  exit  program\nM); 


/*Get  menu  choice*/ 
i=getch(); 
i=tolower( i); 
_clearscreen(_GCLEARSCREEN); 

/,vSelect  memory  option*/ 
if(i='a'){ 
s  e lect_mem( MEM_A ) ; 
goto  loop; ) 

/*Select  memory  option*/ 
else  if(  i='b'  )  { 
select_mem(MEM_B); 
goto  loop; } 

/*Display  memory  option*/ 
else  if(i='c  ){ 
display_mem(MEM_A); 
goto  loop; ) 

/*Display  memory  option*/ 
else  if(i='d^){ 

display_mem(MEM_B); 
goto  loop; ) 

/*Continuous  acquisition  and  display  option*/ 
else  if(i='e')( 
grab(-l); 
goto  loop;  ) 

/*Acquire  and  display  one  frame  option*/ 
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else  if(i='f'){ 
snap(WAIT); 
goto  loop; } 

/,vClear  screen  to  black  option*/ 
else  if(i='g' ){ 
sclear(O); 
goto  loop; } 

/*Clear  screen  to  white  option*/ 
else  if(i='h'  )  { 
sclear(255); 
goto  loop; } 

/*Save   image  to  disk  option*/ 
else   if(i='i'){ 
printf("\n\n   f ilename?\n"); 
scanf("%s" ,name); 
printf ( "comment?\n"); 
scanf ( "%s" , comment); 

if((err=saveim(0,0,511,511,EIGHT_BIT,name,comment))<0){ 
switch(err)  { 

case  FILE_ERROR: 

printf  ("error  opening  file\n"); 
break; 
case  FORMAT_ERROR: 

printf ("unknown  format  selected\nM); 
break; 
default: 

printf ("unknown  error%d\n"); 
break;  } 
exit(l);  } 
goto  loop;  j 

/*Read  image  from  disk  option*/ 
else  if(i='  j'){ 
printf ("\n\n  f ilename?\n"); 
scanf ("%s" j name); 

if ( (err=readim(0,0,511,511,name,comment))<0) { 
switch(err) { 

case  FILE_ERROR: 

printf ("error  opening  file\n"); 
break; 
case  FORMAT_ERROR: 

printf ("unknown  file  format\n"); 
break; 
default: 

printf  ("unknown  error\n"); 
break; } 
exit(l);  } 
printf ("%s\n" , comment); 
goto  loop; } 

/^Change  memory  and  register  base  addresses  in  software  option*/ 

else  if(i='k'){ 
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regbase=OxlOO; membase=OxDOOOOL; 

printf("\n\n  Default  register  base  address  is  0x100.  \n"); 

print f( "Do  you  want  to  change  it?  type  '  y'  or  ,n'\n"); 

nnn=getch( ); 

if(nnn='n' )goto  nextl; 

if(nnn='y' )  [ 

printf("Enter  new  register  base  address,  e.g.  100  \n"); 

scanf("%x",&regbase);  } 

nextl: 

printf("\n\n  Default  memory  base  address  is  OxDOOOOL.  \n"); 

printf("Do  you  want  to  change  it?  type  'y'  or  'n'\nM); 

res=getch( ); 

if(res='n'  )goto  next2; 

if(res='y'){ 

printf ("Enter  new  memory  base  address.   e.g.  D0000\n"); 

scanf ("%lx",&membase);  } 

next2: 

sethdw(regbase,membase,DUAL); 

goto  loop; } 

/*LUT  analysis  option  */ 
else  if(i=U'){ 
_setvideomode( 16); 

select_mem(MEM_A);  display_mem(MEM_A);  grab(-l); 
loopl60: 
_clearscreen(_GCLEARSCREEN); 

/^Display  current  LUT*/ 
lutcoord( ); 

ralut(grp, 0,0, 256, array); 
displaylut( ); 

/"Present  LUT  analysis  options*/ 

_settextposition( 1,1); 

print f("m=mod  lut;  n=mod2;  h=histogram;  e=equalize;  l=linlut; 

i=invlut   "); 
_settextposition(2 ,1); 
printf ("a=abs lut;  p=loglut;  g=grab;  f=freeze;  s=stop;  "); 

/*Get  LUT  analysis  choice*/ 

c=getch(); 

c=tolower(c); 

/^Continuous  acquisition  and  display  option*/ 
if(c=='g' ) (select_mem(MEM_A); 
display_mem(MEM_A);  grab( -1);  goto  loopl60; } 

/*Acquire  and  display  a  single  frame  option*/ 
else  if(c='  f '  )  (select_mem(MEM_A); 
display_mem(MEM_A); snap( WAIT); goto  loopl60; } 

/*Exit  analysis  option*/ 
else  if(c=  s')goto  stop30; 

/*Histogram  option*/ 

else  if(c=='h  )goto  part2; 
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/^Histogram  equalization  option*/ 

else  if(c='e  ) (histequal( ); goto  loopl60;  } 

/*Load  linear  LUT  option*/ 

else  if(c=' 1* )  { linlut(grp,0);  goto  loopl60;  } 

/*Load  inverse  LUT  option*/ 

else  if(c=' i' )  {invlut(grp,0);  goto  loopl60;  } 

/*Load  absolute  value  LUT  option*/ 

else  if(c=' a' )  fabslut(grp,0);  goto  loopl60;  } 

/*Load  logrithmic  LUT  option*/ 

else  if(c='p' )  { loglut(grp,0);  goto  loopl60;  } 

/*Modify  LUT  using  slope  change  option*/ 
else  if (0='^  )goto  mod; 

/*Modify  LUT  by  drawing  option*/ 
else  if(c='m' )goto  parti; 

elsefgoto  loopl60; } 

/*LUT  drawing  routine*/ 

parti: 

pp=180; qq=170;  blx=pp;  elx=blx-l; 

loopl50: 

_settextposition( 1,1); 

printf("b=beginline;  e=endline;  d=drawline;  p=pointer;  s=stop;   "); 

_settextposition(2, 1); 

printfC"  "); 

_moveto(pp-2 ,qq); _lineto(pp+2 ,qq); 

_moveto(pp,qq-2);_lineto(pp,qq+2); 

c=getch(); 

c=tolower(c); 

if(c='3'){pp=pp+10;  if(pp>304)pp=304;  } 


else  if(c='7 

else  if(c='  1 

else  if(c='9 

else  if(c='6 

else  i  f ( c== ' 4 

else  if(c='2 

else  if(c='8 

else  if(c=='b 

else  if(c='e 

else  if(c='s 

else  if(c='d 

else  if(c=='p 

elsefgoto  loopl50; } 

goto  loopl50; 

draw: 

if(blx==elx)elx=blx+l; 

_moveto(blx,bly);_lineto(elx,ely); 

for(  l=blx;  K=elx;  1=1+1)  { 

dely=( float )bly-( float )ely; delx=( float )elx-( float )blx; 


){pp=pp-10;if(pp<51)pp=51;} 
){qq=qq+10; if (qq>299)qq=299; } 
){qq=qq-10; if (qq<46)qq=46;  } 
){pp=pp+l; if(pp>304)pp=304; } 
){pp=pp-l; if(pp<51)pp=51; } 
){qq=qq+l;if(qq>299)qq=299;} 
)fqq=qq-l; if (qq<46)qq=46; } 
){blx=pp;bly=qq; } 
){elx=pp;ely=qq;  } 
)goto  loopl60; 
)goto  draw; 
){pointer();grab(-l);  } 
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f=(300-(float)bly)+((dely/delx)*(l-(float)blx)); 

array[  1-50  =(int)f;  j 
walut(grp,0,0,256,array);  displaylut( ); 
goto  parti; 

/*Slope  changing  routine*/ 

mod: 

_settextposition( 1,1); 

print f("s=stop  "); 

_settextposition(2, 1); 

printf("  "); 

1=0; 

modi: 

if(array[  1+1]  >0){a=l;  goto  mod2;  }  1=1+1;  goto  modi; 

mod2: 

1=255; 

mod3: 

if(array[  1-1]  <255)(b=l;  goto  mod4;  }  1=1-1;  goto  mod3; 

mod4: 

h=a;  g=b;  for(  1=0;  1<=255;  l=l+l)array2[  1]  =array[  1] ; 

mod5: 

c=getch(); 

c=tolower(c); 

if(c=,7')g=g7l; 

else  if(c=' 9' )g=g+l; 

else  if(c=' l' )h=h-l; 

else  if(c='3')h=h+l; 

else  if(c='  s'  )goto  loopl60; 

elsefgoto  loopl60; } 

hf=( f loat)h; af=( f loat)a;  bf=( f loat)b; gf=( f loat)g; 

yi=(255*(hf-af))/(bf-af-gf+hf); 

for(y=0;  y<h;  y=y+l)array[  y]  =0; 

for(y=h;  y<( int)yi;  y=y+l) { 

delta=((hf-af)/(yi-hf))*(yi-(float)y); 

array[ y]  =array2[ ( int )( ( f loat)y-delta)j ; } 
for(y=( int)yi; y<g; y=y+l) { 

delta=((gf-bf)/(yi-gf))*(yi-( float )y); 

array[y]  =array2[  (int)((f  loat)y-delta)j  ;  } 
for(y=g;  y<=255;  y=y+l)array[ y] =255; 
displaylut( ); walut(grp,0,0,256,array); goto  mod5; 

/"Histogram  routine*/ 

part2: 

_settextposition( 1,1); 

printf("d=default;  r=rectangle;  "); 

_settextposition(2,l); 

printf("  "); 

c=getch(); 

c=tolower(c); 

if  ( c='  d' )  { a=40;  b=45;  da=410;  db=405;  } 

else  if(c='  r1  )  {mrectangle(  );  grab(-l);  } 

elsefgoto  part2; j 

linlut( INPUT, 0); linlut(grp,0); 

select_raem(MEM_A); display_mem(MEM_A); snap( -1); 

walut( INPUT, 0,0, 256, array); 

select_mem(MEM_B); display_mem(MEM_B); snap( -1); 
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linlut( INPUT, 0); 

next45: 

_settextposition(  1, 1); 

printf ("l=com  linhist;  m=com  modhist; 

d=dis  linhist;  e=dis  modhist;"); 
_settextposition(2, 1); 

print f("s=stop  "); 

c=getch(); 
c=tolower(c); 
if(c=,l,)comvalsl(); 
else  if(c=='m, )comvals2( ); 
else  if(c='d,)[_clearscreen(_GCLEARSCREEN);  lutcoord(  );  display lut(); 

histcoord();plotvalsl(); linlut( GREEN, 0); display_mem(MEM_A);  } 
else  if(c='e' )  {_clearscreen(_GCLEARSCREEN);  lutcoord(); 

display lut ( ) ; histcoord( ) ; plotvals2( ) ; walut ( GREEN, 0 , 0 , 256 , array) ; 

display_mem(MEM_A);  } 
else  if  ( c= '  s '  )  { linlut( GREEN , 0) ;  walut( grp . 0,0, 256 ,  array) ; 
_settextposition(2, 1); printf ("        ); 
select_mem(MEM_A);  display_mem(MEM_A);  grab( -1); 
goto  loopl60; } 
else^goto  next45; } 
goto  next45; 

/*Exit  analysis*/ 

stop30: 

_settextposition( 1, 1); 

printf ("  "); 

_settextposition(2, 1); 

printfC"  "); 

system("mode  C080"); 

a=30; b=30;  da=450;  db=450; 

goto  loop; ) 


/*Save  LUT  to  disk  option*/ 
else  if(  i='m'  )  { 
printf  ("enter  LUT  file  name\n"); 
scanf ( "%s" , lutnarae); 

save_lut(lutname,grp,0,0,256,NEW_FILE); 
goto  loop; } 

/*Read  LUT  from  disk  option*/ 
else  if  (  i='n'  )  { 
printf ("enter  LUT  file  narae\nM); 
scanf ( "%s" , lutname2); 
read_lut( lutnarae2); 
goto  loop; j 

/*Select  output  LUT  for  loading  option*/ 
else  if  (  i='  o'  )  { 
grp=GREEN; 
goto  loop;  j 

/^Select  input  LUT  for  loading  option*/ 
else  if  (  i='p'  )  { 
grp=INPUT; 
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goto  loop; } 

/*Zoom  option*/ 

else   if(i='q'){ 

select_mem(MEM_A); display_mem(MEM_A); 

printf("\n\n  To  move  upper   left   corner  of  rectangle  type:\n\n"); 

printf('"8'  for  up\n"): 

print f("' 2'  for  down\nf'); 

printf("'4'  for  left\n"); 

print f  ( " '  6 '    for   right\n\n\n") : 

printf("type    'z'    to  zoom\n\n*  ); 

loop70: 

while(!kbhit()){snap(WAIT);  rectangle(a,b,240,235 ,0);  } 

e=getch(); 

e=tolower(e); 

if(e='8'){b=b-10;if(b<0)b=479;  } 

else  if(e='2,)[b=b+10;if(b>479)b=0;  } 

else  if(e='4'  ){a=a-10;  if(a<0)a=512;  } 

else  if(e=V){a=a+10;  if (a>512)a=0;  } 

else  if(e— '  z' )goto  stop3; 

elsefgoto  loop70; } 

goto  loop70; 

stop3: 

pan(a); scroll(b); zoom(MAG_2); 

spot20: 

print f("\n\n  type  'm'  to  return  to  menu\n"); 

c=getch(); 

c=tolower(c); 

if(c=V  ){da=240;db=235;goto  loop;  } 

else(goto  spot20; }   } 

/,vUnzoora  option*/ 

else  if  (  i='  r1  )  ( 
zoom(MAG_l);  scroll( -b); pan( -a); 

select_mera(MEM_B);display_mem(MEM_B); sclear(255); 
select_mem(MEM_A); display_raera(MEM_A); sclear(255); 
initregs( ); set lut(grp,0); 

select_mem(MEM_A); display_mem(MEM_A); 

grab(-l); 

a=40; b=45; da=410;  db=405; 

goto  loop; ) 

/^Sharpen  option*/ 

else   if( i==' s' ) { 

print  f(  "\n\n\n\n\n\n\n\n\n\n\n 
PROCESSING   IMAGE\n"); 
sharpen( a , b , da , db , POS ITI VE ) ; 
goto   loop; } 


/*Low  pass   filter  option*/ 
else   if(i='t')( 
print  f  ( "\n\n\n\n\n\n\n\n\n\n\n 
PROCESSING   IMAGE\nM); 
lopass ( a, b, da ,db, POSITIVE); 
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goto   loop; } 


/*Sobel   edge   filter  option*/ 
else   if( i=='u' ) { 
print  f  ( "\n\n\n\n\n\n\n\n\n\n\n 
PROCESSING  IMAGE\n"); 
sobel(a,b,da,db, 16); 
goto   loop; } 

/*Image  average  option*/ 
else  if(i=V){ 
print f( "number  of  frames ?\n"); 
scanf("%d",&num); 
_clearscreen(_GCLEARSCREEN); 
print  f  ( "\n\n\n\n\n\n\n\n\n\n\n 
PROCESSING   IMAGE\n"); 
average( 10, 10,490,470,num,num); 
goto   loop; j 

/*Framegrabber  off  option*/ 
else  if(i=V) 
fgoff(); 

} 

/*This  function  places  a  moveable  pointer  on  video  screen  and  feeds 
position  and  intensity  level  to  computer  monitor*/ 
pointer( ) { 

select_mem(MEM_A);  display_mem(MEM_A); 

_settextposition( 1,1); 

prmtfC    s=stop  ); 

_settextposition(2, 1); 

printfC  ); 

ppl=255;qql=240; 

loop50: 

snap(WAIT); 

intens=rpixel(ppl,qql); 

_settextposition( 2 , 65); 

printf  ( "x=%d\n"  ,ppl); 

_settextposition( 3,65); 

printf  ("y=%d\n"  ,qql); 

_settextposition(4,65): 

pr int  f ( "  int ens  ity=%d\n'  ,  intens  ) ; 

if ( intens<=127)color=255; 

if ( intens>127)color=0; 

xl=ppl-4; x2=ppl+4; yyl=qql; y2=qql; x3=ppl; x4=ppl;  y3=qql-4;  y4=qql+4; 

clr_line(xl,yyl,x2,y2,color); clr_line(x3,y3,x4,y4,color); 

while(!kbhit())jsnap(WAIT);  clr_line(xl,yyl,x2,y2, color); 

clr_line(x3,y3,x4,y4, color); } 

c=getch(); 

c=tolower(c); 

if(c='3')[ppl=ppl+10;  if(ppl>511)ppl=0;  } 

else  if(c=^7'  )fppl=ppl-10;  if (ppl<0)ppl=511;  } 

else  if(c='  l'  )  (qql=qql+10;  if (qql>479)qql=0;  ) 

else  if(c='9'  )  {qql=qql-10;  if  (qql<0)qql=479;  } 
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) [ppl=ppl+l; if (ppl>511)ppl=0 
){ppl=ppl-l;if(ppl<0)ppl=511 
) jqql=qql+l; if (qql>479)qql=0 
) [qql=qql-l; if (qql<0)qql=479 
)goto  loop50; 
)goto  stop; 


else  if(c='6 
else  if(c='4 
else  if(c='2 
else  if(c=*8 
else  if(c='5 
else  if(c='s 
else{goto  loop50; } 
goto  loopSO; 
stop: ; 

_settextposition(2,65); 
printfC"  fi); 

_settextposit ion( 3 , 65) ; 
printf("  f'); 

_settextposit ion( 4 , 65 ) ; 
prlntf("  ");} 

/*This  function  places  a  rectangle  whose  size  and  location  may  be 
changed  on  the  video  screen.   Size  and  location  parameters  are  then 
used  in  other  functions*/ 
mrectangle( ) { 

select_mera(MEM_A);  display_mem(MEM_A); 

_settextposition( 1,1); 

printf ("s=stop  "); 

_settextposition(2, 1); 

printf(  ); 

loop60: 

while(!kbhit()){snap(WAIT); rectangle(a,b,da,db,0);  } 

e=getch(); 

e=tolower(e); 

if(e='8')fb=b-10;  if (b<0)b=479;  } 

else  if(e='2'  ){b=b+10;  if (b>479)b=0;  ) 

else  if(e='4'  ){a=a-10;  if(a<0)a=512;  ) 

else  if(e=V  ){a=a+10;  if (a>512)a=0;  j 

else  if(e='9'){da=da+10;  if(da+a>512)da=512-a;  } 

else  if(e=l7')[da=da-10;  if (da<10)da=10;  ) 

else  if(e='3'){db=db+10;if(db+b>479)db=479-b;  ) 

else  if(e='l'){db=db-10;if(db<10)db=10;  ) 

else  if ( e==' s ' )goto  stop2; 

else{goto  loop60; j 

goto  loop60; 

stop2:  ;  } 

/"This  function  performs  histogram  equalization  and  allows  for  manual  or 
automatic  updates*/ 

histequal( ) { 
loop3: 

_settextposition( 1, 1); 

print f(Md=default;  r=rectangle;  "); 

_settextposition(2, 1); 

printf("  "); 

c=getch(); 
c=tolower( c); 

if  ( c='  d' )  { a=40;  b=45;  da=410;  db=405;  } 
else  if(c='  r' )  {mrectangle( );  grab(  -1);  j 
else{goto  loop3; } 
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loopll: 

_settextposition( 1, 1); 

print f("c=continuous  update;  m=manual  update;  "); 

_settextposition(2, 1); 

printf ("  "); 

c=getch(); 
c=tolower(c); 
if(c='c' )goto  loopl2; 
else  if(c=lm,); 
else{goto  loopll; } 
color=4; 
loop2: 

select_mem(MEM_B);  display_mem(MEM_B); 
snap(WAIT); 
select_mera(MEM_A); 
display_mem(MEM_A); 
grab(-l); 

s  e lect_mem( MEM_B ) ; 
histeq(grp,0,a,b,da,db); 
ralut(grp,0,0 ,256, array); 
_setcolor(color  %   16); 
display lut( ); 

color=color+5;  if (color>9)color=4; 
loop4: 

_settextposition( 1,1); 
printf("u=update;  c=change  param;  f=freezeframe;  g=grab;  s=stop;"); 
_settextposition(2, 1); 

printfC  ); 

c=getch(); 
c=tolower(c); 
if(c='u' )goto  loop2; 

else  if(c=' c' )  {_setcolor(9  %  16);  goto  loop3;  } 
else  if(c='  f '  )  (snap(WAIT);  goto  loop4;  } 
else  if(c==' g1 ) |grab( -1); goto  loop4; } 
else  if(c==' s1 ) {_setcolor(9  %  16); goto  stop25; ) 
else^goto  loop4; ) 
loopl2: 
color=4; 

_settextposition( 1, 1); 

printf ("enter  time  between  updates  (integer  2-9)  "); 

_settextposition(2, 1); 
printf (M  "); 

C=B6t ch( ) ' 

if(c='2')vt=2;else  if(c=' 3' )vt=3;  else  if (c='4'  )vt=4; 

else  if(c='  5  '  )vt=5;  else  if(c=' 6' )vt=6; 

else  if(c='7')vt=7;else  if  (c==' 81  )vt=8;  else  if(c='9'  )vt=9; 

elsefgoto  loopl2; ] 

_settextposition( 1, 1); 

print f("c=change  param;  s=stop;  "); 

_settextposition(2, 1); 

printf("  "); 

while(!kbhit()){ 

time(&ltime); 

xx=ltime; 

select_mem(MEM_B); display_raem(MEM_B); 
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snap(WAIT); 

s e 1 ect_mem( MEM_A ) ; 

display_mem(MEM_A); 

grab(-l); 

select_mem(MEM_B); 

his teq( grp , 0 , a , b , da , db) ; 

ralut(grp, 0,0, 256, array); 

_setcolor( color  %  16); 

display lut( ); 

color=color+5;  if(color>9)color=4; 

while(( ltime-xx)<vt)time(&ltime); } 
c=getch(); 
c=tolower(c); 

if(c='c')[_setcolor(9  %  16);  goto  loop3;  } 
else  if(c='  s'  )  (_setcolor(9  %   16);  goto  stop25;  } 
else{_setcolor(9  %   16); goto  loop3; } 
stop25: 

se lect_mem( MEM_A) ; 
display  mem(MEM_A); 
} 

/*This  function  displays  the  current  LUT  shape  on  the  computer  monitor*/ 
display lut() { 

c=_getcolor( ); 

_setcolor(2  %   16); 

_moveto( 50 , 300) ;  _lineto( 305 , 45 ) ; 

_setcolor(c  %  16); 

_moveto(50,300); 

for(l=50; 1<=304;  1=1+1) ( 

_lineto( 1+1, 300-array[ 1-49] ); }  } 

/*This  function  draws  the  coordinates  for  the  histogram  display*/ 
histcoord( ) { 

_moveto(370,300);_lineto(370,40);_moveto(370,300);_lineto(630,300); 

_settextposition( 23,47); printf("%d" ,0); 

_settextposition(23,58);printf("%d",100); 

_settextposition(23,71);printf(M%d",200); 

_settextposition( 15 ,43); print f ("%d" ,2); 

_settextposition(8,43); printf("%d",4); 

for( 1=370; 1<=620; 1=1+50) { 

_moveto( 1, 300 );_lineto( 1,305); 
_moveto(370,670-l);_lineto(365, 670-1);  }  ) 

/*This  function  computes  a  histogram  for  the  contents  of  memory  A*/ 
comvalsK )  { 

_settextposition( 1,1); 

printf("  ^  "); 

_settextposition(2, 1);  printf("        "); 

_sett extpos it ion( 13,50); print f ( "COMPUTING  HISTOGRAM" ) ; 

select_mem(MEM_A); 

histogram(a,b,da,db, 1, l,0,histvalsl); 

_settextposition( 13,50); printf("  "); } 

/*This  function  computes  a  histogram  for  the  contents  of  memory  B*/ 
comvals2( ) { 

_settextposition( 1,1); 
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printf(  ); 

_settextposition(2,l);printf(  ); 

_settextposition( 13,50);  printf ("COMPUTING  HISTOGRAM"); 

s e lect_mem( MEM_B ) ; 

histogram(a,b,da,db, 1, l,0,histvals2); 

_settextposition( 13,50);  printf ("  ");  } 

/*This  function  plots  the  cotnvalsl  output  histogram*/ 
plotvalsl(){ 

_moveto(370,300); 
for( 1=370;  1<=625; 1=1+1) { 
_moveto( 1,300); 
_lineto(l,300-(histvalsl[ 1-370] /20)); }} 

/*This  function  plots  the  comvals2  output  histogram*/ 
plotvals2()[ 

_moveto(370,300); 
for( 1=370;  1<=625; 1=1+1) { 
_moveto( 1,300); 
_lineto(l,300-(histvals2[ 1-370] /20)); }} 

/*This  function  draws  the  coordinates  for  the  LUT  display*/ 
lutcoord( ) { 

_moveto(50,300);_lineto(50,40):_moveto(50,300);_lineto(310,300); 

_settextposition(23,7); printf(  '%d",0); 

_settextposition(23,18);printf("%d",100); 

_settextposition(23,31);printf("%d",200); 

_settextposition( 15,3); printf (  %d" , 100); 

_settextposition(8,3);printf("%d",200); 

for( 1=50; 1<=300; 1=1+50) { 

_moveto( l,300);_lineto( 1,305); 
_moveto(50,350-l);_lineto(45, 350-1); }] 
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